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Abstract. Since its discovery by Max Planck in 1900, the Planck constant h has been 
demonstrated to be an universal constant, and its numerical value has been accurately 
determined based on experiments. Up to the present however the physical origin of this 
fundamental constant has not been well understood, and the numerical value of it has not 
been ab initio predicted, h is characteristic in two respects: 1) it is a universal constant 
with respect to all (quasi-) stationary dynamical processes of all matter particles and 
radiation fields, and 2) it has a specific numerical value. A theoretical inference of h, 
and a corresponding accounting for the physical origin of h, therefore needs be achieved 
in both respects. This paper presents a theoretical exploration in the first respect, a 
mathematical inference of the universal constancy of h, based on the second law of 
thermodynamics, the principle of least action and the probability theory. 



1. Introduction 

Max Planck hypothesised energy quantum in 1900 [1] in order to resolve the then 
demonstrated large discrepancy between classical kinetic theory prediction and experimental 
data for black body radiation in the high frequency end, known as the "ultraviolet 
catastrophe" . That work laid the foundation of quantum mechanics. Through the formal 
quantum mechanics subsequently developed by E Schrodinger, W Heisenberg, P Dirac and 
others in 1920s-30s along with a wide variety of experiments subsequently performed, it has 
become well established that h is a universal constant for all (quasi-) stationary processes 
of all matter particles and radiation fields. 

At the scale h, a classically point matter particle turns to behave as a wave, W. At 
this scale, the motion of a particle i of mass Mj and potential energy V is, in directly non- 
relativistic terms here, governed by the the Schrodinger equation, with a,kn = t**, a^n = h, 

ia kn d t V = HV, H = -^V 2 + V. (1) 

Electromagnetic radiation fields may be subjected to the same equation (1), provided one 
regards, as is commonly done in quantum electrodynamics, (1) as describing the motion of 
a charged harmonic oscillator emitting the radiation. In a similar way, (1) presents also a 
governing equation for the generation of the total internal radiation fields, and hence the 
total energy and mass of a particle (see Sec. 3). Our central concern in this paper is the 
fundamental Planck constant h. Equation (1) and other basic quantum-mechanical laws, on 
the other hand, present the context within which h is meaningfully defined. It will suffice 
that we refer to (1) as representing the overall quantum-mechanical laws, since from (1) 
or the solutions to (1), the basic dynamical variables of a particle and their relations such 
as the thermal and the total eigen energies and momenta, and the Heisenberg uncertainty 
relations, among others, can all be derived. 
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It is well understood that the wave equation (1) in essence describes the stationary state 
of a particle of a wave function ^(x, t) and probability density p(x) = \&(x, t)\ 2 . This follows 
directly from the fact that every equation (1) is associated with a continuity equation for p, 
jfi + V(pv) = 0, where v is the velocity of p. In a confined space, the dynamical variable 
solutions to (1) for the particle are furthermore generally quantised, with respect to the 
Planck constant h as prescribed in (1). Evidently, however, a quantised stationary-state 
solution will result irrespective of what the akn value is. akn can in fact be an arbitrary 
real-valued parameter for each specified energy level n and (positional) degree of freedom 
k. More generally, akn can be dependent on the nature of the particle system and its 
environment, as Planck similarly had suggested in [lb]. 

Today, h is taken for granted to be a universal constant. Its numerical value has 
hitherto been determined by fitting such theoretical (or empirical) functions as the blackbody 
radiation spectrum, and photo-electric work function, etc., to experimental data. These 
appear to have all resulted in (a/s„ =)h — 6.6260Q x 1CP 34 Js. The physical origin of the 
h in respect both to its universal constancy and its specific value, however, has hitherto 
remained an open question. There is no known logical reason why all afc„'s need be equal, 
in particular equal to h. 

The universal constancy of h is also responsible for a basic concept in statistical 
mechanics, namely, in the phase space of, say, TV particles having D positional degrees 
of freedom each, there exists a smallest volume, h DN , which can be accessed by each 
microscopic state, or microstate, of the TV-particle system. This concept was also originally 
introduced by Planck [lb], based on an hypothesis of absolute entropy and a presumption 
of the validity of the postulate of equal a priori probabilities, namely that all microstates 
of an isolated thermodynamic system in equilibrium are equally accessible. Similarly, there 
is hitherto no known logical reason why the involved postulate, which is mainly in question 
here, holds except that based on it, statistical-mechanical solutions for thermodynamic 
systems have shown overall agreement with experiment. 

An understanding of the origin of the universal constancy of h will shed light on one 
of the key issues (another being the origin of energy quantisation) regarding the origin of 
quantum mechanics; and this is also a simultaneous step needed for an ab initio prediction 
of the numerical value of h. In this paper, starting with a thermodynamic system with its 
component particles described by the wave equation (1) each, with akn assumed an arbitrary 
parameter, we shall infer the postulate of equal a priori probabilities (Sees. 2-3) and the 
universal constancy of akn, hence h (Sec. 4), based on the second law of thermodynamics and 
Maupertuis-Jacobi's principle of least action, severing as our first principles here, combined 
with the probability theory. 

2. Maximum entropy 

We consider, until Sec. 3, iV (identical) particles in an enclosure of volume V, that are (a) 
isolated (from the environment), (b) weakly interacting only, and (c) distinguishable. Here, 
(b) is said in the sense that the particles may exchange energies but present no correlations in 
respect to occupancy of states. The TV-particle system is in (or approaching) thermodynamic 
equilibrium and is thus, owing to condition (a), of (or approaching) a fixed internal energy 
U. The particles are assumed to be each governed by the wave equations (1), with a^ n 
assumed arbitrary parameters. We shall regard the particles as intrinsically waves here; at 
the large quantum number (n) limit these will reduce to the classical "point" particles. 

Assume that the TV particles have each D (positional) degrees of freedom, or dimensions. 
The quantum state of a particle i (i = 1,2, . . . , N) at time t is thus completely specified 
by D eigen wave functions <Pfc n (xk,i)'s, with k — Di — (D — 1), Di — (D — 2), . . . , Di. In 
a confined space [0, Lk] (e.g. Lk — V 1 ^ 3 for a free particle in V = Lf), the eigen energy 
Ekn of the particle i is in general quantised, with n = 1,2,..., 9Tfc , where n is the principal 
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quantum numberj and Vtko is the total number of possible energy levels in fcth dimension. 
There are a total Y^k=Di-(D-i) number of distinct possible states for particle i. 

The microstate of the TV-particles at a given time t is accordingly completely specified 
by D x TV cigen wave functions of the TV individual particles. Each ith particle at time 
t statistically lies in a definite energy level n, of an energy E kn w, with n suffixed for 
cxplicitness here by the particle index (i). And the TV-particles has a definite total internal 
energy U = Y^iLi J2k=Di-(D-i) ^fcnW = Z)i=i E mM ■ Of the total Vlko levels for ith particle 
in fcth dimension, there can be 9^ (< Vtko) levels that satisfy condition (a), given by n's 
lying in a narrow range [n a , rib], defining therefore an TV-particle system with fixed N,V, U. 
The corresponding possible microstates may be expressed by the set of eigen wave functions 

» = {#l„(i) » • ■ • ! ^Dn(i) I • • • ! *fcn<*> . • • ■ 5 ^(fllV-(fl-l))n(») . • • • . %W)n(*0 < « W < } 

= {Vl v ,*2 v , ■ ,%W)> = 1, 2, . . . ,ft} (2) 

which make up a microcanonical ensemble. In the second of Eqs. (2), each microstate 
specified by the original set of level indexes nW,...,nW,...,nW is re -labelled, in an 
arbitrary order relative to the original, by a single running index v, v being an unique value 
for each distinct microstate. ft denotes the total number of possible microstates satisfying 
the condition (a), and is formally given as 

N Di N Di 

n = Yl 9tfc = 5Z ft *' Qi= 2 mk - ^ 

i=l k=Di-{D-l) i=l k=Di-(D-l) 

We assume here, as condition (d), that fi >> TV so that the situation of more than one 
particle occupying the same state will be rare. The ft, or Ot^ wave functions of a particle 
i span a Hilbert space in the state- or energy- representation, where the state vector 

IS &k,ens = Y^(% = i c knW^knW, or &i,ens = Hk=Di-(D-l) &k,ens = Z^(0=l c inW&inW i 

=Di-(D-i) ^fc,n< s > ■ The ft wave functions (2) span a Hilbert space Jf? of 

the TV particles, where the state vector is W ens = YliLi ^i,ens — X^=i c ^^i &v = 
(Sign)y / r!f r i iy (xi) • ... • &( DN ) v (x DN ) being the total wave function of the TV independent 
and (for r = 1) distinguishable particles here. The c v etc. in the above are the amplitudes 
of states v, etc. 

The microscopic states of the TV-particles may be alternatively described in an usual 
2DTV-dimensional phase space & spanned by DN space coordinates x\, . . . , Xdn and DN 
momentum coordinates Pi,.. . ,Pdn of the TV particles. A vth distinct microstate is in 
completely specified by a volume element (Aari„, . . . , Axmm v ; AP\ V , . . . , AP/ DN \ V ) located 
about a fixed point (x\ v , . . . , X( DN \ V \ P\ v , . . . , PmN)v), which are bounded between adjacent 
energy levels each as will be formally defined by (23)-(24); and this has a one to one 
correspondence with a point B„, = {tfq„, . . . , &(dn)v}i m ^ ■ Accordingly, each z/th volume 
element occupies a volume 

DN DN 

K = \\ Ax kv AP kv = JJ a kv , a kv = Ax kv AP kv , v = 1, 2, . . . , Q; (4) 

fe=l k=l 

and all of the Ct volume elements occupy a total volume B = X^=i b v in 3? ■ By its definition 
(4), b v is the smallest volume accessible to a state v. 

% This holds also for a particle described by a set of quantum numbers, e.g. n, I, m, s of a bound atomic 
electron. Here, I (= 0, 1, . . . , n — 1) and m = ±1 reflect the geometric orientations of the total angular 
momentum (nh) and total magnetic moment of the orbiting electron in magnetic field, s describes the spin 
which is a permanent property of the electron and, as of an any elementary particle, is never subjected to 
exchange. The l,m,s each yield no independent or direct contributions to entropy, and thus need not be 
the concern of this paper. 
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At any instant of time t the N individual particles statistically lie in their respective 
definite eigen states <F ln (i) , . . . , ^(£>jvW w ) wrtn tne probabilities |c ln (i> | 2 , • • • , \cmN) n w ? '■ 
And the iV-particles lie statistically in a definite microstate M v , or volume b v in with 
the probability |c„| 2 . Over long time, as the result of particle-particle interactions under 
condition (a), the TV- particles will statistically explore all the accessible possible microstates. 
To observe all of the ft possible statistical states — providing also accessible — at least once 
requires a measurement to be made over long time. To facilitate a measurement required 
only to be made at any one time, we may instead employ a Gibbs ensemble consisting of 
E(<ft) replicas of the original iV-particle system, 

G = {&i lt ,&2v l ,...,*kn,---,V(DN) li \n = 1,2,...,-} (5) 

which all have different microstates and yet identical macroscopic properties, the fixed 
N,V,U here. 

We now want to accommodate the E replicas of the iV-particle system in the phase space 
J 2 *, by imagining these as E objects "thrown" into 8?, which will spontaneously attain their 
equilibrium positions after a relaxation time. Consistent with our presumption that the 
afc„'s are not necessarily all equal (Sec. 1), however, we do not a priori assume here, as the 
postulate of equal a priori probabilities instead does, that all the ft possible microstates B 
are equally accessible to the iV-particles. Or equivalently, we do not assume that all the ft 
volume elements in & can be readily uniformly occupied by the E replicas, and this leads 
to E 7^ ft (or E < ft). The actual accessibility or occupancy of B by G, hence the actual 
correspondence between IB and G, is to be determined. 

The degree of accessibility of a vth state is directly proportional to the probability 
y v for the iV-particles to be found in the ;/th state upon a measurement at time t, or 
equivalently, for a replica \i to be found in the i/th volume element in 3 21 . We assume that 
the 2£W-dimensional phase space in respect to a replica pertains to a geometric nature 
in the same sense as the position space in respect to, say, a dart in the dart game. Thus 
to attempt to accommodate or "throw" a replica \i in a volume element v in region B in 
is like to attempt to throw a dart into a small surface element of area <j(r,dr,d9) on 
a dartboard D. Based on dart game experiment (see e.g. [2]) performed in an effectively 
homogeneous space (assuming gravity field etc. is negligible), the larger the area a is, the 
more likely the dart will hit <i on 2). The probability T a for a successful hit at the target 
area a is thus y a oc a. Making direct analogy, granted with a homogeneous phase space in 
region B where U is everywhere the same, then the larger the volume b v of the vth volume 
element is, the more likely a replica \x will "hit", or be accommodated in, b v . That is, the 
probability 3>„ for a successful accommodation of a replica fi in b v is oc b v . Multiplying b v 
by 1/B gives the normalised probability 

9 V = ^, i/ = l,2,...,Q. (6) 

Evidently, ? u - (1/B) £? =1 b v = 1. 

Now according to the second law of thermodynamics, the disorder liness, hence the 
entropy S, of an isolated system in (or approaching) thermodynamic equilibrium is (or 
approaches) maximum. Evidently, the disorderliness of the TV-particles will be maximum 
if the system explores as uniformly as possible all of the ft possible microstates over time. 
This thus requires: 

All the ft possible microstates are at any time t simultaneously equally accessible to the 
iV-particle system. (7) 

(7) may be restated in terms of the E replicas and the phase space as 

The ft volume elements in are simultaneously uniformly occupied by 
with one and only one replica occupying one volume element at a time. 



S = ft replicas, 

(8) 
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The simultaneously accessible states in (7), or the simultaneous occupancies of the 
possible states stated less abstractly in (8), represent simultaneous events that are, due to 
condition (b), independent with one another. The independence here removes the possibility 
of correlated occupancies of B by the replicas due to particle interactions. The condition (d) 
restrains two or more replicas /U, fi' from occupying the same state v at the same time. And 
the ft value is specified by the number of eigen solutions of (1) for the N particles under the 
fixed N,V, U condition. On the other hand, it remains possible that some regions in region 
B in 8? are densely populated and some other regions rarely as a result of the geometric 
property of the phase space, so that effectively S < ft; this situation dissatisfies (8) or (7). 
If such situation is to be avoided, the phase space must possess certain quality as may be 
identified as follows: The probability for the occurrence of the f> simultaneous independent 
events (7) or (8) is according to probability theory given by 

(7) or (8) will be maximally fulfilled if 7 

ens is maxirnuin; and ^P ens is a maximum (T ens . max ) 

if 

sy ens = o, 5 2 y ens < o. (io) 

To (universally) determine y e ns.max by solving (10) using calculus will inevitably be 
unfeasible, since the volume elements in & can be in arbitrary shapes owing to the diverse 
shapes and dynamics of physical systems and Tens is therefore not generally expressible 
into a universal analytical function. A maximum y ens solution however can be obtained 
readily by means of algebraic method as follows. The T ens given by (9) is equivalent to the 
geometric mean of CP^'s 

^=(U^ = (n "=^ }1/ " , (ii) 

V— 1 

raised to the power of ft, i.e. T ens = . On the other hand, the arithmetic average is 

It follows from the theorem of inequality of the arithmetic-geometric means that T v < 

(3V) = (jy), and hence T ens =T V < (3\/) ft (= ($y) a ), ^ tne ^V s > being nonnegative real 
functions, are not all equal. And, 

^ u .max — {^f) — ^ens.maxi — ^u.max ) ' 

i.e. IP„ and 7 en s are maxima each, if and only if the nonnegative real functions JVs are all 
equal, !Pi = ?2 = ■ • • = CPj>- Or equivalently, (14) holds if and only if, on substituting (6) 
for 3>„, all b„s are equal to one another and hence to a common value denoted by b , 

h = b 2 = ■ ■ ■ = ba = b - (15) 

With (15), (4b) becomes B = Y^=i = ^b and (6) becomes 
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(16) states that, for the TV-particles under conditions (a)-(d), all the ft possible microstates 
are equally accessible, which is just the statement of the fundamental postulate of equal a 
priori probabilities. The entropy of the N- particles follows to be a maximum, given as 

S = -k B ^2y u \ny„ = k B \nil^2- = k B \nil. (17) 

I/=l v=\ n 

The maximum probability CP„ and entropy S of (16)-(17) are achieved each in a statistical 
sense, as are the statements (9) and (10), and hence allow for statistical fluctuations. 

The formal steps from (9) to (14) in the above have yielded a few significant equations 
(11), (12) and (14). Without these we could have as well obtained (15) and (16) by continuing 
the geometrical argument as made for (6). Namely, because T u oc b v , all TVs are thus equal, 
whence satisfying (7) or (8), if all 6„'s are equal. 

3. Global system 

For defining Ct and b u , and for arriving at the basic conclusions in Sec. 2, Eq. (16) in 
particular, the N particles need not be all identical. The N particles may be of different 
species, chemical compositions and/or energy forms (see below), and they may be of different 
interaction and motion schemes. These may for example include iV a atoms, N e electrons 
(free or bound), N p ht photons (electromagnetic radiation energy quanta), N p h n phonons 
(sound energy quantum) , and so forth. Although, the total quantities of the mixed particle 
system need be given by the sums over all particle species, etc., e.g. N — J2 a= iN a , 
ft = J2 a N a D a yi a . And mixed particles will typically involve different degrees of freedom 
(DOF's). If two different DOF's represent the same energy, such as the energy levels of 
a harmonic oscillator and the quanta of its radiation fields, evidently only one should be 
considered for entropy. If a DOF is internal, such as the orbiting motion of an atomic 
electron, this DOF should not be (directly) considered for entropy. 

In usual applications, Eq. (1) describes the kinetic motion of a matter particle such as 
an electron in an applied potential field V, or indirectly the thermal radiation fields. The 
Ei n , Pin etc. are the dynamical variables associated with such motions. Alternatively, (1) 
also describes the motion of a charge, of a dynamical mass 5DT g j (in the place of Mi), in 
a (quadratic) vacuum potential field across a distance b v <~ 10~ 18 m, which generates in 
terms of the IED model [4] the total, internal motion of a simple matter particle such as 
an electron, proton (see a systematic treatment in [4]c). (E in -^)E qin — nhHi = nMiC 2 is 
then the total energy and Pj„ the total linear momentum of the particle, with fii = ^jr-, 
and Mj (with n = 1) the mass of the resulting stable IED particle. The basic conclusions of 
Sec. 2 apply evidently directly to this latter system — a particle of an alternative (total as 
contrasted to thermal) energy form. 

Equation (16) has been otherwise arrived at for a system of N particles under the 
conditions (a)-(d) that are also the conditions for the direct application of the postulate 
of equal a priori probabilities in statistical mechanics. So (16) may be extended to other 
(important) variant systems following the same well-established "correction" procedures (see 
e.g. [3]) in statistical mechanics, as are briefly outlined in (i)-(iv) below, (i) Two important 
variations from the "isolated system" of cond. (a) are the "closed isothermal system", 
whose U may fluctuate and N, V, T are fixed, and " open isothermal system" , whose U, N 
may fluctuate and V,T,n (where ft is the chemical potential) are fixed. The U and U,N 
fluctuations yield inhomogeneities in their 3Vs are now given by the canonical (or 
Boltzmann) and grand canonical distributions 



y{U) =Q{U)e~ 1 ^ T /Q(T,V,N), T(U, N) — Q(U, N)e w s r ~ /Z(V, T, fi), (18) 



7 



where Q = r J2u> ft(U')e V , and Z = rJ2u>,N< ^( u> > N ') e " > with r = 1 and 1 / Nl 
for N distinguishable and indistinguishable particles. Separately, the iV-particle system may 
change energy, by AC/, by interacting with another system by means of thermal interaction 
under fixed external parameters, or by adiabatic change of external parameter(s), or by a 
combination of the two. In any of the cases, ft is unchanged, (ii) The N particles may 
be correlated with one another, hence a deviation from cond. (b). The usual method to 
extend beyond (b) is virial expansion of the equation of state, (iii) The N particles may 
be indistinguishable with one another, hence a deviation from cond. (c). This is especially 
the feature with many quantum particles confined in a common small geometry, where their 
wave functions mutually overlap to a high degree. The standard correction to (c) is, as is 
similarly done in (18), by means of dividing out Nl number of indistinguishable states from 
ft, giving fti.D = ■ T ma Y be very low such that ft ^ N and the situation of more 
than one particle occupying the same states becomes frequent, hence a deviation from cond. 
(d). The standard approach to such system is the replacement of the Boltzmann statistics 
by quantum statistics (Fcrmi-Dirac or Bose-Einstein statistics) where restrictions are made 
to the number of accessible states by applying symmetry properties of the wave functions of 
particles. Of the systems of (i)-(iv), the b v and ft values are unchanged. Only the extents 
or forms of the contributions of their ft to T v may be modified and thus corrected for. The 
corrections may be appreciated as being to scale the systems back to obeying conditions 
(a)-(d). A recognition of this feature permits one to generalise the approach of "scaling 
back" also to other possible variant systems not yet considered. 

The diverse forms of particle systems of this section and Sec. 2 together resemble 
a major important part of the "global system" of the physical world in which energy 
transformations and exchanges in units of integer times an energy quantum hx angular- 
frequency (e.g. the fli earlier) have been observed. To this global system the inference of 
the postulate of equal a priori probabilities in Sec. 2, Eq. (16), is valid in the fashion as 
case by case discussed above. 

4. Least action 

Insofar as all the &„'s are ensured equal (to bo) and thus S(U) is maximum for a prescribed 
ft value, it is irrelevant that what the a kv values, and hence the Ax kv and APt„ values 
according to (4), are unless otherwise constrained. In the following we shall not assume any 
pre-established knowledge of the a ku 's (esp. their being equal to h), but shall investigate 
the characteristics of the afc„'s, or the akn's below, regarding their constancy based on 
Maupertuis-Jacobi's least action principle 

Of the ./V particles of Sec. 2, we focus now on an individual particle (i) in motion 
at a component velocity v k in fcth dimension confined in a space interval [0, Lfc]. The 
particle is assumed (quasi-) stationary and described by an eigen wave function and 
probability density \^ kn {x k , t)\ 2 governed by Eq. (1). The so-described particle is generally 
extensive across [0, L k ] at any time t, and oscillating over all time at any location x k in 
[0, L k \. Its dynamical variables, such as linear momentum P kn (x k ,t) — Miv' kn , kinetic 
energy X' kn (x k ,t) = \P kn {x k ,t)v' kn {x k ,t), potential energy V kn (x k (t)), and Hamiltonian 
E 'kn( x k-,~t) = 3C' fe „(a;fc,i) + V kn (x k (t)), are accordingly distributed functions in [0, L k ]. 

For this undulatory extensive stationary particle, there exists in general a characteristic 
space interval AX kn , such that only over AX kn , or [0, AX kn ], the particle can be wholly 
and meaningfully defined, or in fact exists in stationary state when condition for complete 
constructive/destructive interference presents. For a variety of monochromatic periodic 
particle processes, AX kn = nAXji; and the minimum of AX kn , AX k \, is given by the 
distance of one-full cycle of the periodic process in space, which is e.g. one wavelength of 
a plane de Broglie wave (this however does not apply to a plane electromagnetic wave 
of one energy quantum, or a photon, which is not directly governed by Eq 1), or 4x 
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the amplitude of a harmonic oscillation. There thus exists a corresponding characteristic 
time interval AT kn = AX kn /(v' kn ) needed for the particle to traverse [0,AX kn ]. And in 
the case of AX kn = nAXji, AT kn — nAT k \\ the minimum of AT kn , AT k \, is equal to 
one full period of time. AX kn for a particle at a given energy level n must in turn be 
accommodated in an external environment, such as between enclosure "walls" at x k = 
and L k . Hence AX kn = AX kn (L k ). For example, for a free particle of a wave function 

■ / 2-tt E n j.\ 

\P n (x,t) — Ce A d n x »n ' and wavelength Adn confined in a one-dimensional box [0,L] 
(with the subscripts k omitted here), we find AX n = nAXi, AXi = Adi- But Ad n = ^ as 
given by the boundary conditions tf'(O) = ^{L). Hence AX n (= nAd n ) = 2L. 

We define for the particle a Maupertuis-Jacobi's action integral§ in [0, AX kn ] as, 

A kn = PUx k ,t)dx k = , kn{ k \ ' dx k = / 2X' kn {x k ,t)dt, (19) 

JO JO v kn\ X k: t ) Jo 

where dx k /v' kn — dt. (19) may be re-written as 

Akn = PknAX kn , A kn — (2X kn )AT kn , (20) 

where P kn ,0C kn are the expectation values of P kn ,X' kn or of the corresponding operators 
P kn ,X kn acting on \P kn (with W kn assumed normalised in [0, AX kn ]): 

^ i>AX kn rAXkn 

Pkn= A y / P'kn( x k,t)dx k =1 >P kn (x k ,t)P kn <P kn (x k ,t)dx k , 

LA^kn Jo Jo 

3C kn = -t-=- / X' kn (x k ,t)dt = / W* kn {x kl t)X kn ^ kn {x k ,t)dt- E kn = X kn +V kn .(21) 

LAJ-kn Jo Jo 

(Pkn and X kn of the typical stationary particle systems in applications, and E kn of any such 
systems, are in fact independent of x kl t.) We further define the difference (Maupertuis- 
Jacobi's) action integral as the absolute difference between A kn 's of two adjacent levels 
n + 1 and n: 

rAx kn 

akn = \A kn+ i - A kn \; or a kn = AP kn Ax kn = / AP kn (x k ,t)dx k , 

Jo 

r-At kn 

a kn = (2AX kn )At kn = / 2AX' kn (x k ,t)dt = AE kn At kn , (22) 
Jo 

where Ax kn and At kn are the mean characteristic space and time intervals defined by 

a _ \Pkn+\AXkn+l — PknAXkn\ * , _ \^kn+l ATkn+l 

AXkn = T~~r> ' M kn = T^? 5 ^ 23 > 

APkn AX kn 

and 

AP kn = \P k ( n +i) - Pkn\, AX kn = \X kn+ i — X kn \, AE kn — \E k ( n +i) — E kn \; (24) 
AX kn = AV kn , thus AE kn = AXkn + AV kn = 2AX kn . (25) 

AP kn and AX' kn are similarly defined. As it may be checked against the typical quantum 
systems in applications, using Ax kn ,At kn defined in (23) as combined with the eigen 
solutions for P ni X n will indeed yield the exact forms of Heisenberg relations. For a spatially 
confined particle, the AP kn , AX kn , a kn , etc, are quantised as are the P kn ,X kn and E kn 
(Sec. 2) and A kn - Moreover, the dk n , AP kn , AX kn , etc are in the above each positively 
defined since their magnitudes only, not their signs, will be relevant. The mechanical energy 

§ A kn may be re-written as A kn = f[X' kn - V fc ' n + %' kn + V^Jdt = f[£' kn + E' kn ]dt, where C' kn = X' kn - 

is the Lagrangian function. J £' krl dt = S kn defines a distinct action integral which is widely in use today 

and which, as may be easily seen, is yet unsuited for the present problem. 
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difference AE kn corresponds to an energy quantum of harmonic radiation field, or a photon, 
emitted upon transition from n + 1 to n, for which the relations (25) are always valid. 

It is by the requirement of the least action principle that the A kn and accordingly a kn 
need be minimum each and hence satisfy the conditions 

(a) : 5A kn = 0, S 2 A kn > 0; (b) : Sa kn = SA kn+1 - SA kn = 0, S 2 a kn > 0. (26) 

Since Ax kn , At kn , AP kn and AE kn are finite for the AX kn and AT kn being finite — 
intrinsically for an intrinsically extensive IED particle[4], (26) must have nontrivial, finite 
valued solutions. Since Y^n'=o akn ' = •^■fcm to achieve a minimum A kn solution for (26a) 
it suffices, ultimately for the final solution form, that each a kn is minimised according to 
(26b). 

For N particles in the generalised sense of the global system of Sec. 3, least action 
should ultimately be satisfied by all of the a kn a) 's of all energy levels of all N particles, and 
hence by the sum of a kv 's over all v values: 

DN DN DN DN 

I v = J2 akv = Yl ^ p kvAx kv = 1AX kv At kv ; 5I V = 5^a kv =Q 1 5 2 I U > 0. (27) 

k=l k=l k=l fc=l 

Similarly as for 7 V earlier, a general differentiable function for J„ for all particles in the 
global system, which in general may be of arbitrary dynamics and enclosure geometries, 
and thus a general solution for (27b) given in terms of calculus, will be absent. However, a 
minimum I v can be readily obtained using algebraic method as follows. 

In (4) we have already written down a product equation of all a k „'s, with the a kv 's 
positively defined. With (15), (4) is rewritten as 

DN 

K = Y[a k v = b , u = 0,1,..., St. (28) 
fe=i 

For a chosen global system of a total fixed N number of particles in or approaching 
equilibrium with fixed macroscopic properties, the product of the DN parameters a kv 's, 
bo, given in (28) is a fixed value. Then, according to the theorem (a corollary of the theorem 
of inequality of the arithmetic-geometric means) for finding extrema, the sum I v of the DN 
positive a kv 's as given by (27a) is a minimum if and only if all the positive a kv 's are equal 
to one another, and hence to a common value denoted by a , 

DN 

a\ v = a<i v = . . . = d( DN)v = ao, v = 1,2, ...,ft; I v . m in — ^ Qq = DNa . (29) 

k 

Substituting (22b)-(c) in (29a)-(b) gives 
AP fel Ax fc i = . . . = AP kv Ax kv = AE kl At kl = ... = AE kv At kv = a , v = 1, 2, . . . , ft; 

DN DN 

AP kv Ax ku = Y 2AX kv At ku = DNa . (30) 

k=l k=l 

Substituting (29) in (22a) gives (for A n+ i > A n ) A kn+ \ = A kn + a . Calculating ,Afc ra +i 
successively, we obtain A k \ = A k o + a , A k 2 = A k \ + a = A k o + 2a , and so on; and 

A kn (= P kn AX kn = 2K kn AT kn ) = A ka + na (31) 

Conversely, if the sum /„ of the positive a kll 's is set to a fixed value, concretely the minimum 
Iv.min, then according to the foregoing theorem stated in reverse order, the product of the 
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positive cifcj/s, b (= ^k=i a kv), is a maximum if and only if all the positive a^'s are equal, 
as given in (29). With (29) and (30), (28) is written as 

DN DN DN 

bo = ]J ak " = II Ap k»&x ku = Y[ 2AX kv At kv = a° N = b v . max . (32) 

k=\ k=l k=l 

Except for being required as general and as inclusive as possible, the global system may 
be chosen differently at different occasions or times, whose fe ' s will thus be different. As 
long as the different global systems share at least one common particle — a trivial condition 
which can be readily satisfied in both theory and practice, however, the same universal ao 
must be preserved for all choices of the global systems at all times, ao thus is universal 
irrespective of the choice of the global system and time. And the foregoing inference of the 
universal constancy of ao is valid to the extend that the second law of thermodynamics, to 
which no violation has hitherto been observed, and the least action principle are valid. The 
latter, the least action principle, too is a general dynamical law as much as the former, in the 
sense that from it, such basic mechanical laws as the Euler-Lagrange equations, Newton's 
second law, the Schrodinger equation, and the de Broglie relations (see e.g. [5]), can each 
be derived. 

By its having the dimensions "joule x second", its universal constancy as 
mathematically inferred in the foregoing, and its basic relationships with particle dynamical 
variables as given e.g. by (1), and (30)-(32), a is therefore identifiable with the Planck 
constant h, a — h. Accordingly (1) is identifiable with the Schrodinger equation, (30a) 
with the Heisenberg uncertainty relations 

AP kl Ax kl = ... = AP kv Ax kv = AE kl At kl = ... = AE kv At kv = h, (33) 

(31) with e.g. the de Broglie relations P kn A kn = nh, E kn r kn = nh (with Ajo = 0, 
AX kn = A kn , AT kn = r kn , and n effectively continuous) of a free particle, and (32) with the 
usual equation for the smallest volume h DN accessible to a statistical-mechanical microstate 
in phase space. 
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Sept, 2012, where the author has enjoyed interesting conversations with several participants, 
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radiation quantum) at the workshop, to P-I Johansson for private financial support of the 
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